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Abstract 



A new tight-binding total energy method, which has been shown to accurately 
predict ground state properties of transition and noble metals, is applied to 
Manganese, the element with the most complex ground state structure among 
the d metals. We show that the tight-binding method correctly predicts the 



■ ground state structure of Mn, and offers some insight into the magnetic prop- 



erties of this state. 
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I 

>■ . Most elements in the periodic table crystallize in the fee, bec, hep and diamond structures. 

Among the few exceptions is Manganese, which has an equilibrium structure, denoted a-Mn, 
which contains 29 atoms in the unit cell 0. First principles total energy methods, such as 
the full-potential Linearized Augmented Plane Wave (LAPW) method |2||| are not very 
efficient in such systems, especially since the a-Mn phase has five internal parameters which 
must be adjusted to minimize the total energy in order to calculate the correct structure at 
each volume. Calculation of elastic constants, phonon frequencies, surface energies, vacancy 
formation energies, and other properties require even more computational effort. A reliable 
approximate method, based on first-principles results, is necessary for efficient computational 
study of complicated crystals such as Manganese. 

In a recent paper, Sigalas and Papaconstantopoulos Q introduced the idea that the en- 
ergy bands of Augmented Plane Wave (APW) calculations for cubic structures at different 
volumes can be fit to a non-orthogonal tight-binding (TB) Hamiltonian whose matrix ele- 
ments are functions of the distance between pairs of atoms. The sum of eigenvalues resulting 
from the above TB Hamiltonian, together with a pair potential, were used to fit the total 
energies of the APW calculation, thus obtaining an interpolation formula that was employed 
to calculate the total energy for non-cubic structures. This procedure was applied to cal- 
culate the elastic constants of Pd, Ir, Au, Rh and Ta, which showed fairly good agreement 
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with the experimental values. The phonon spectra and density of states for Au were also 
calculated [f|, again in reasonable agreement with experiment. 

In a subsequent paper, Cohen, Mehl and Papaconstantopoulos || made dramatic im- 
provements to the above approach. They eliminated the pair potential in the fitting of 
the total energy; employed environment-dependent on-site TB parameters; and introduced 
exponentially damped polynomial expansions of the hopping and overlap integrals, thus ex- 
tending the parametrization to an arbitrary number of neighbors. This new total energy 
methodology |J was applied to calculate elastic constants, phonon spectra and vacancy 
formation energies for the noble metals and other transition metals. The results were im- 
pressive. Starting from only fee and bec structures, the method correctly predicted the 
ground state structure in all of the elements tested, including those which exhibit an hexag- 
onal close-packed (hep) ground state. 

This paper shows how the new TB method can be applied to Manganese. We first 
performed paramagnetic LAPW calculations at five volumes in each of the monatomic fee 
and bec structures, and then determined a set of TB parameters which reproduced 
the electronic structure and total energy of these structures. We then used the resulting 
Hamiltonian to compute the total energy of Manganese in the a-Mn, /3-Mn f cc > bec, 
hep, and simple cubic (sc) structures. The a-Mn structure has a bec unit cell containing 
atoms, with five internal parameters. Its space group is 743m-T| [0]. The twenty- nine atoms 
are divided up into four types, with all atoms of a given type equivalent by symmetry. There 
is one atom of type J, located at the origin; four atoms of type II, located at a(x±, xi, X\) and 
equivalent points, where a is the cubic lattice constant; twelve atoms of type III, located at 
a(x2, X2, z-i) and equivalent points, and twelve atoms of type IV, located at a(xs, X3, z 3 ) and 
equivalent points. The /3-Mn structure has a simple cubic structure containing twenty atoms, 
with two internal parameters. The space group is P4 3 32-0 6 KJ. There are eight atoms of 
type /, located at a(xx, x 1: xi) and equivalent sites; and twelve atoms of type 77, located 
at o(l/8, x 2 , 1/4 + x 2 ) and equivalent sites. Neither structure has an inversion site, so we 
must solve the generalized eigenvalue problem for Hermitian matrices. The large number of 
atoms, coupled with the necessity of minimizing the total energy with respect to the internal 
parameters, makes the determination of structural properties difficult to handle by first- 
principles total energy electronic structure calculations. Within the tight-binding method, 
however, the calculation is relatively easy. In Table | we show the equilibrium volume, 
relative energy, and bulk modulus as calculated by our TB procedure for Manganese, as well 
as Technetium and Copper for comparison. Fig. [I] shows the energy/volume relationship 
for several of the lower energy phases, and Table compares the equilibrium structural 
parameters with the experimental ones. Our TB Hamiltonian correctly predicts the ground 
state structure of Manganese. The calculation shows that the /3-Mn phase is close in energy 
to the a-Mn phase, indicating that it is a likely candidate for the high-temperature phase 
of Manganese, in agreement with experiment J1J. Our calculations also predict that a-Mn 
will transform into an hep structure at a pressure of 50 GPa. 

Of course this agreement with experiment could be an artifact of the way we constructed 
the Hamiltonian. To test this, we constructed TB Hamiltonians for Technetium, which is also 
a column VIIB element, and Copper by determining Tight-Binding parameters which repro- 
duced the results of APW (for Technetium) and LAPW (for Copper) calculations as outlined 
above. The resulting equation of state data is presented in Table |I[ Our parametrization 
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correctly predicts that the hep phase is the ground state for Technetium although the a-Mn 
and /3-Mn phase are close in energy. Copper is correctly predicted to be in the fee phase, 
while the a-Mn and (3-Mn phases are respectively 6.0 mRy and 5.4 mRy higher than the fee 
phase. 

Our calculations for Manganese were performed assuming a paramagnetic phase, while 
experiment || and theory fl- [l~l|1 suggest that most phases of Manganese exhibit some form 
of magnetism, and that the a phase is antiferromagnetic [HPfl. This does not affect the 
fact that our method demonstrates that the a phase is the ground state of Manganese, 
because the addition of magnetism can only lower the energy of the a phase. Magnetism 
will, however, affect the volume of the ground state phase. Our equilibrium lattice constant 
is about 6% smaller than the experimental lattice constant. This error can be partially 
attributed to the neglect of magnetism and partially to the error inherent in the LDA . 
Our calculated internal parameters (Table |J) are almost identical to the experimentally 
measured parameters, so we conclude that the internal parameters are not changed by 
magnetism. 

Since the current formalism is not set up for spin-polarized calculations, we use our 
paramagnetic TB Hamiltonian to calculate the electronic density of states (DOS) for the 
a-Mn phase. Fig. |^ shows the total DOS of a-Mn as well as the d partial DOS for the four 
different atom types. The width of the d states appears the same for all sites. However, 
there are differences in the details of the DOS structure. In particular, the Fermi level 
values of the DOS differ substantially. This is shown in Table |TXT] where we note that the 
first two sites have DOS values which are a factor of two larger than the other two sites. 
Also from Table [III] we note that the s and p-like DOS at Ep are very small. We then 
applied the Stoner criterion [13,14] using the DOS at the Fermi level in conjunction with a 
matrix element derived from fee and bec calculations [TTS|. Using an approximate value of 
If = 0.03Ry []15 |, we obtained values of the Stoner nl of about 0.7 for atoms on sites / and 



77, and about 0.4 for atoms on sites /// and IV. This is consistent with first-principles 
band-structure calculations for the moments on the atoms ||, where it is found that atoms 
I and // have large moments, but atoms /// and IV have smaller moments. 

We conclude that our TB total energy method is capable of predicting the correct total 
energy ordering of various structures, including the complicated a-Mn structure, with com- 
putational costs orders of magnitude lower than standard first-principles calculations. In 
addition this scheme provides reliable energy bands and DOS for all phases of Manganese. 

We wish to thank R. Cohen, W. Pickett, and D. Singh for useful discussions. 
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FIGURES 



FIG. 1. The total energy of Manganese in several structures as a function of volume, obtained 
from the tight-binding Hamiltonian outlined in the text. We show energy per atom versus volume 
per atom for ease in comparison. The o symbols indicate the fee and bec phase LAPW energies 
used in the fit. The a-Mn phase is correctly predicted to be the ground state. 

FIG. 2. The electronic density of states (DOS) of the a phase of Manganese, as well as the d 
density of states for the four distinct sites. The dotted vertical line represents the Fermi energy. 



5 



TABLES 



TABLE I. The equilibrium volume (Vo) in Bohr 3 and energy (Eq) in mRy per atom, and bulk 
modulus (-Bo) in GPa for Manganese, Technetium, and Copper, as calculated by the Tight-Binding 
method. Eq is set to zero for the ground state energy of each element. 



Atom 




fee 


bec 


hep 


sc 


diamond 


a-Mn 


/3-Mn 


Mn 


v 


68.7 


68.7 


68.4 


74.0 


96.8 


69.3 


69.6 




E 


8.2 


15.5 


3.1 


91.0 


172.0 


0.0 


1.6 




B 


315 


324 


314 


199 


119 


320 


318 


Tc 


Vo 


94.2 


95.4 


93.6 


102.5 


127.2 


95.1 


95.1 




Eq 


6.5 


23.9 


0.0 


57.2 


73.2 


0.2 


2.6 




Bo 


309 


306 


303 


244 


175 


299 


298 


Cu 


Vo 


73.6 


73.9 


73.8 


82.6 


106.5 


75.2 


75.5 




E 


0.0 


3.5 


1.5 


25.1 


70.9 


6.0 


5.4 




Bo 


190 


186 


186 


141 


64 


176 


177 



TABLE II. The experimental and TB equilibrium lattice and internal parameters for the a 
structure of Manganese. The bec unit cell has twenty-nine atoms, divided into four classes as 
explained in the text. 

a(A) Xl X2 Z2 %3 23 

Experiment 8.9129 0.31765 0.35711 0.03470 0.08968 0.28211 

Tight-Binding 8.41 0.31719 0.35787 0.03964 0.08971 0.27983 



TABLE III. The total electronic density of states at the Fermi level for a-Mn at the minimum 
energy volume predicted by the tight-binding calculations. The partial DOS for the s, p, and 
d states are shown for each atom type (see text). The Stoner criterion parameter is calculated 
assuming Ip = 0.03 Ry [15]. The coordination numbers are those assigned by Donohue [0. 



Type 



Atom 

Number Coordination 



DOS at E F 
Total Partial 
(States/Ry/Unit Cell) (States/Ry/Atom) 

s p d 



Stoner 
Criterion 



I 

II 
III 

IV 



1 

4 
12 
12 



12 
10 
13 
11 



458.39706 0.07778 0.70738 22.81580 0.68 

0.27634 0.43595 25.07444 0.75 

0.21981 0.64905 14.41710 0.43 

0.13672 0.60299 11.61058 0.35 
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DC 
LU 



0.000 



■0.005 - 



■0.010 - 



■0.015 - 



■0.020 - 



■0.025 - 



■0.030 - 



60 




65 



70 
V (au 3 ) 



75 



80 



Mn 



ENERGY (eV) 



1400.0 



700.0 - 




0.1 0.3 0.5 

ENERGY (Ry) 



